Contact Dependence of Carrier Injection in Carbon Nanotubes: An Ab Initio Study 
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We combine ab initio density functional theory with transport calculations to provide a micro- 
scopic basis for distinguishing between good and poor metal contacts to nanotubes. Comparing Ti 
and Pd as examples of different contact metals, we trace back the observed superiority of Pd to the 
nature of the metal-nanotube hybridization. Based on large scale Landauer transport calculations, 
we suggest that the 'optimum' metal-nanotube contact combines a weak hybridization with a large 
contact length between the metal and the nanotube. 
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A major challenge linked to the use of carbon nan- 
otubes |l| in future electronic devices is to understand 
the profound effect of the nanotube-metal contact on 
transport. Weak nanotube-metal coupling, found in nan- 
otubes deposited on metal electrodes, has been shown to 
cause Coulomb blockade behavior |2j. In spite of sig- 
nificant progress in maximizing the contact area by de- 
positing metal on top of nanotubes Q , the transparency 
of such contacts exhibits strong sample-to-sample vari- 
ations and depends strongly on the contact metal. Re- 
ports of low contact resistance between nanotubes and 
Au or Au/Cr 0,0 are in stark contrast to the high re- 
sistance observed in nanotube contacts with Au/Ti 
The transparency of Pd-based contacts has been reported 
as superior in comparison to using Ti, Pt and Al as 
contact metals 0, H, @- Additional modulation of the 
Pd-nanotube contact transparency has been reportedly 
achieved by modulating the gate voltage [lfl |. Reports 
suggesting that carrier injection occurs only at the edge 
of the contact region Q appear to contradict the observed 
dependence of the contact resistance on the length of the 
contact @. 

Published theoretical results include studies of the elec- 
tronic structure at a nanotube- Au interface and trans- 
port properties of a nanotube- Al junction [J]J. Ab initio 
calculations furthermore suggest that Ti contacts may 
be superior to those with Al or Au 12], and that the 
Schottky barrier between semiconducting tubes and Pd 
is lower than with Au or Pt 13]. Due to the limitation 
to specific contact geometries and small system dimen- 
sions, however, general trends are hard to extract, and 
an extrapolation to experimentally relevant system sizes 
is difficult. 

Here we combine ab initio electronic structure stud- 
ies with large scale transport calculations to gain micro- 
scopic insight into the relative importance of the inter- 
face morphology, the type of the contact metal, and the 
length of the contact region when optimizing the metal- 
nanotube contact. Ab initio density functional studies 
were used to determine the charge redistribution and 
electrostatic potential in the contact region. In a second 



step, the electronic structure results were mapped onto 
a model tight-binding Hamiltonian suitable for transport 
calculations. We found that transmission is maximized 
in the case of weak metal-nanotube coupling, exhibited 
by extended Pd contacts. 

To gain insight into the electronic structure in the 
contact region, we performed density functional theory 
(DFT) calculations of Ti and Pd monolayers interacting 
with a graphene layer. We described valence electrons by 
Troullier-Martins pseudopotentials and used the Perdew- 
Zunger form of the exchange-correlation functional in the 
local density approximation (LDA) to DFT, as imple- 
mented in the SIESTA code [l4|. With a double-zeta 
basis and a 100 Ry energy cutoff in the plane-wave ex- 
pansions of the electron density and potential, we found 
the total energy to be converged to <1 meV/atom. We 
performed a full structure optimization to determine the 
equilibrium adsorption geometry, the adsorption energy, 
and the local charge redistribution caused by the metal- 
graphene interaction. Since the interatomic distances in 
bulk Pd (2.7 A) and Ti (2.95 A) lie close to the honey- 




FIG. 1: (Color online) Charge density redistribution Ap(r) = 
pMe/c( r ) — PMe(r) — pc(r) in (a) Pd and (b) Ti monolayers 
interacting with a graphene layer, indicating regions of charge 
depletion and excess with respect to the superposition of iso- 
lated layers, (c) Schematic double-layer geometry in top view, 
with the cutting plane used in (a) and (b) indicated by the 
dash-dotted line. 
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comb spacing in graphene (2.46 A), we only considered 
epitaxial adsorption. For both Pd and Ti, we found a 
slight preference for the sixfold hollow site on graphite, 
depicted schematically in Fig.^c). For Pd, we found the 
equilibrium inter-layer distance to be 3.2 A, consistent 
with a relatively weak, mostly covalent bond energy of 
0.3 eV per Pd atom. The interaction between an epitax- 
ial Ti monolayer and graphene was only insignificantly 
stronger with 0.4 eV per Ti atom at an inter-layer dis- 
tance of 3.0 A. 

The quality of nanotube-electrode contacts has been 
shown to depend sensitively on the Schottky-barrier in 
semiconducting nanotubes [l5| and band bending in 
metallic nanotubes, both reflecting the charge transfer 
within the junction. Our Mulliken population analysis 
indicates a net charge transfer of only psO.I electrons from 
Pd and Ti to the graphene layer. More useful informa- 
tion is contained in the charge redistribution, depicted 
in Fig.Q] Results for Pd electrodes, shown in Fig. ^a), 
suggest an accumulation of excess charge in the region 
between Pd and graphene layers. As seen in Fig. db), 
the charge redistribution in Ti/graphcnc is very differ- 
ent, suggesting charge accumulation in the atomic layers 
and depopulation of the interlayer region and thus an 
increase of the interlayer scattering potential. The lower 
scattering potential and the populated interlayer state at 
the Pd/C junction appear well suited for carrier injection 
into the nanotube, making the Pd/C contact superior to 
the Ti/C contact. 

Besides the charge transfer, the contact quality de- 
pends even more sensitively on the nanotube-metal hy- 
bridization [J^, which is addressed in Fig. [21 for the 
Pd/graphene system. The density of states at E-p as- 
sumes a large value, which is a pre-requisite for a good 
contact. The electronic band structure of the system, de- 
picted in Fig. |2ta) , suggests that all states are closely re- 
lated to either Pd or graphene states. In the Pd/graphene 
system, the graphene bands are rigidly shifted by AEc = 
0.374 eV and the metal bands by A£: Pd = -0.020 eV 




FIG. 2: (Color online) (a) Electronic band structure E(k) of 
a Pd monolayer interacting with a graphene layer, (b) Details 
of E(k) in (a) near the Fermi level, denned as Ep = 0. 



with respect to the isolated layers. In the Ti/graphene 
system, the rigid band shift at the junction is much 
stronger, AE C = —1.15 eV, and has the opposite sign 
to Pd. 

Especially interesting for the transparency of the con- 
tact is the nature of Pd-C rehybridization, which is best 
visible in Fig. 13b) close to the Fermi level. Particu- 
larly clear is the hybridization between Pdc? z 2 and Cp z 
orbitals, which causes a «0.15 eV band splitting about 
0.5 eV below £f, in the vicinity of the K point. Since 
K denotes also the Fermi momentum of graphene, this 
occurrence of Pd-C hybridization near this fc-point sug- 
gests an efficient way to inject carriers into graphene near 
the Fermi energy without involving phonons to conserve 
momentum. 

To obtain quantitative information about the effect of 
the junction geometry and hybridization on the trans- 
parency of the contact, we performed large scale quan- 
tum transport calculations of nanotubes in contact with 
metal electrodes within the Landauer approach . Our 
calculations for nanotube segments exceeding 10 2 nm 
were facilitated using an efficient O(N) decimation al- 
gorithm (lif . Our transport calculations were based on a 
simplified tight-binding Hamiltonian, describing only the 
interaction between Pdd z 2 and Cp z orbitals. We found 
that the electronic band structure of the Pd/graphene 
system near Ep, depicted in Fig.|2£b), can be reproduced 
by considering the hybridization between the ppir band of 
graphene, associated with 70 = 2.66 eV, and a Pd-based 
band, using tpd/c ~ 0.15 eV for the hopping integral 
between Pd and each of the six C neighbors. Such a 
simple mapping turned out insufficient to describe the 
hybridization between Ti and graphene. Based on typ- 
ical band repulsion observed in that system, the Ti-C 
hopping integral ixi/c ~ 0-3 eV should be about twice 
that found for Pd. 

Results of our transport calculations for metal- 
nanotube junctions are summarized in Fig. [3J In our 
schematic view of an extended contact, depicted in 
Fig. Eta), we distinguish three regions within a finite 
tube. The central region of length L$, describing an un- 
perturbed nanotube, is connected at both ends to contact 
regions of varying length L c . In both extended contact re- 
gions, each atom is coupled to the coating metal electrode 
in a similar way, as previously considered in Ref. |lf)j . 

In the model examined in the following, we chose di- 
agonal wide-band leads, which contacted each atom of 
the nanotube independently. In this case, the coupling 
is described by an energy- independent, purely imaginary 
self energy T,(E) — — iA, where A = ij^gA/Me^F) can 
be extracted from our ah initio results. Using the calcu- 
lated density of states Mmc(Ef) at Pd and Ti surfaces, 
which is of the order 1 eV -1 , we obtain Apd ~ 0.06 eV 
and Axi ~ 0.3 eV, assuming that each carbon atom is in 
direct contact with three metal atoms. 

Transmission T through a molecular conductor is gen- 
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FIG. 3: (Color online) (a) Schematic geometry of the (6, 6) 
nanotube in contact with metal leads, used in the calculation 
of the contact reflection coefficient p. (b) Contact reflection 
coefficient p as a function of the nanotube-metal coupling A 
and the contact length L c . (c) Cuts through the contour plot 
(b) at selected values of A, showing p as a function of L c . 
The effective contact length is emphasized by a heavy 
solid line in (b) and by data points in (c). 



erally limited by the number of channels iV c h, which de- 
pends on the band structure of a perfectly periodic sys- 
tem as T < Tbands = -^ch- The reason for the effective 
transmission through a conductor with contacts T(L C , A) 
being lower than through the ideal infinite system Tbands 
lies in the reflection at the contacts. To quantify the qual- 
ity of a nanotube contact, we define the contact reflection 
coefficient by 



p(Tc,A) 



1 



1 



(T(L C ,A)> Tba 



(1) 



The average is taken around the Fermi level, in a re- 
gion between the first van-Hove singularities, similar to 
a transmission convolution capturing hot electron ef- 
fects [2(| ■ Our results in Fig. |3| based on ensemble aver- 
aging in Eq. agree quantitatively with those obtained 
by averaging the transmission using the Fermi distribu- 
tion, for a wide range of temperatures. Physically, the 
total resistance R of such an idealized system could be 
separated as R — -Rband + R c , where i?band = l/(2Go) 
is given by the quantum limit of two open channels with 
conductance Go = 2e 2 /h each. Assuming zero temper- 



ature and neglecting disorder effects, R c ~ p/Gq origi- 
nates only from sub-optimal contacts and may theoreti- 
cally be arbitrarily small. 

We generally expect zero transmission in the limiting 
cases of vanishing contact length, L c = 0, and vanishing 
coupling, A = 0. For finite values of L c and A, how- 
ever, it is not obvious if a combination of strong coupling 
and short contact is superior to a combination of weak 
coupling and a long contact. To obtain a quantitative 
answer, we calculated p for a contact to a (6, 6) armchair 
nanotube with Lq = 100 nm as a function of L c and A. 
We found p to be independent of the tube diameter, as 
long as the energy range used in the averaging avoids 
subband-related van Hove singularities. We also found 
no dependence on Lq, as long as Lq was much larger 
than a unit cell size. Consequently, p should be only a 
function of L c and A. 

Our results for p(L c , A) are depicted in Figs. 3(b) and 
(c). We find the contour plot of p in Fig.^Jb) separated 
into two regions by a line of "minimum contact reflec- 
tion" L° ff (A). fn short contacts, the transparency is re- 
stricted by L c and p increases with decreasing L c due to 
a generalized Breit-Wigner broadening of the resonances. 
The pronounced ripples found at small values of A and 
L c are not numerical artefacts, but rather reflect the in- 
terplay between resonances in a finite nanotube segment 
and conduction electrons propagating with the Fermi mo- 
mentum. For very weak coupling, we find {T) E oc T C A. 
For larger values of T C A, however, dissipation along the 
contact region modifies this simple behavior, yielding 
p ex exp (— aL c A). 

In sufficiently long contacts, defined by L c > L\ , p 
becomes independent of L c and is given by p — PqA 2 . 
We found that our results can be reproduced well by 
using po = 0.016 eV~ 2 . The independence of p from 
L c in long contacts is seen clearly in the plots of p(L c ), 
depicted in Fig.^c) for selected values of A. Particularly 
intriguing appears our result that reflection is minimized 
in case of a weaker specific coupling, provided the contact 
is sufficiently long. 

This physical origin of this unexpected behavior is 
schematically illustrated in the insets in Fig. |3c) . Sat- 
uration transparency is reached for relatively short con- 
tacts in the case of strong coupling. In this case, however, 
the abrupt change in the electronic structure between the 
uncoated and the coated nanotube segment causes extra 
reflection. This change is less abrupt in the case of weaker 
coupling and a long effective contact region T° , reducing 
the overall reflection. 

To estimate the effective contact length L° s , we make 
use of the above described expressions for p in the adja- 
cent regions in the (L c , A) plane. The line, where these 
two functions intersect, corresponds to the line of mini- 
mum p, and is given by the analytical expression 



Lf (A) =4 



■ A A 



(2) 
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Here, £ uc — 2.46 A is the unit cell length. The parameters 
a-y = 1.34 eV and «2 = 9.14 eV were obtained by fitting 
our numerical data. In the specific case of Pd and Ti, we 
found Lf (A Pd/c ) w 30 nm and Lf (A Ti/c ) w 4nm. 

Since realistic metal-nanotube contacts are rarely epi- 
taxial, we have considered various forms of disorder in the 
contact region and determined their effect on transport. 
We modelled weak to moderate disorder by randomly 
perturbing the metal-nanotube coupling A with respect 
to AMe/o We did not find noticeable change in p even 
for perturbations of A as large as its reference value. As 
another extreme case, we modelled strongly diluted con- 
tacts by randomly suppressing the interaction between 
individual metal and nanotube atoms in the contact re- 
gion, down to one percent of contacts with respect to the 
epitaxial case. Finally, we considered various forms of 
non-diagonal contributions to the self-energy, modelling 
the cross-coupling between neighboring metal atoms. In 
all the studies, which addressed deviations from epitaxy 
at the interface and our description of the leads, we found 
the same behavior as depicted in Fig. and analytically 
described by Eq. with possibly modified numerical 
values of po, a,\ and ui |18| . 

Our main conclusion, which proved to be robust with 
respect to variations in the details, is that each con- 
tact can be characterized by an effective contact length 
Lf, which depends only on the local metal-nanotube 
coupling, not on the diameter of the tube. Assuming 
that the effective contact length between the nanotube 
and the electrode exceeds L° , which is likely the case in 
most experimental setups, then a higher contact trans- 
parency is expected when the metal-nanotube coupling 
is weak [l^. Especially in very long contacts, the sensi- 
tive p = poA 2 oc ^Mc/c dependence of p on the coupling 
strength £mc/c niay be taken as an important guideline, 
suggesting to minimize coupling and maximize contact 
length to achieve a high contact transparency. In the 
specific case of Pd and Ti contacts, the weaker nanotube- 
metal coupling in the case of Pd is a good explanation 
for the superiority of Pd-based nanotube contacts. 

In conclusion, we combined ab initio density functional 
theory with transport calculations to distinguish micro- 
scopically between 'good' and 'poor' metal contacts to 
nanotubes. Comparing Pd and Ti as examples of dif- 
ferent contact metals, we traced back the observed su- 
periority of Pd to the nature of the metal-nanotube hy- 
bridization. Based on large scale Landauer transport cal- 
culations, we suggest that the 'optimum' metal-nanotube 
contact generally combines a weak hybridization with a 
large contact length of typically few hundred nanometers 
between the metal and the nanotube. 

Of immediate interest is, of course, the general valid- 
ity of our results. We plan additional studies address- 
ing the effect of non-epitaxial contacts and the nature 
of charge carriers |l8| . Particularly interesting in this 
respect should be studying spin injection from ferromag- 



netic contacts and Andreev reflection at the contact to 
superconducting electrodes. 
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